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We propose a double exchange model to describe the RuC>2 planes of RuSr2(Eu,Gd)Cu208- The 
Ru +5 ions are described by localized spins, and additional electrons provided by the superconducting 
CUO2 planes are coupled ferromagnetically to them by Hund rules coupling. We calculate the 
spin structure factor, magnetic susceptibility and magnetization as a function of magnetic field and 
temperature, using a Monte Carlo algorithm in which the Ru +5 spins are treated as classical. Several 
experiments which seemed in contradiction with one another are explained by the theory. 



PACS numbers: 75.30.Vn, 75.10.-b, 75.40. Mg 

RuSr2(Eu,Gd)Cu208 is a very interesting material be- 
cause of the coexistence of magnetic order and super- 
conductivity below Ts ~ 45 KBcl The system orders 
magnetically at Tm = 133 K, but the type of magnetic 
order remains controversial. The first measurements in 
Gd samples indicated ferromagnetic (FM) ordering of the 
Ru momentsjj but neutron diffraction experiments found 
superlattice reflections consistent with an antiferromag- 
netic (AF) order withjjiearest-neighbor spins antiparallel 
in all three directions.El The Ru contribution to the mag- 
netic susceptibility at temperature T > Tm can be ver 
well described by x = C/(T - 6), with 6 = 100 ± 3 K 
The fact that > seems difficult to reconcile with AF 
order at low temperatures. 

From previous intensive research on the similar sys- 
tem YBa2Cu306+:r|j it seems clear that each of the su- 
perconducting Cu02 planes has doping slightly less than 
0.1 holes per Cu ion in order to lead to the observed su- 
perconducting critical temperature T<j. Since Eu or Gd 
are expected to be in the oxidation state +3, assuming 
as usual that the inactive SrO layers are neutral, charge 
balance imposes that the Ru02 planes have an electron 
doping n < 0.2 per Ru ion.B Neglecting covalency with 
O atoms, this means that there are 1 — n Ru +5 and 
n Ru +4 ions per unit cell. The superconducting Cu02 
planes and the magnetic Ru02 planes can be regarded 
as separate entities as a first approximation. This is sup- 
ported by symmetry considerations and band-structure 
calculations We propose that the electronic structure 
of the Ru02 planes can be described by an effective dou- 
ble exchange model, with one localized spin at each site 
representing the Ru +5 ions, and additional n itinerant 
electrons per Ru coupled ferromagnetically with the lo- 
calized spins. 

The model is the two-dimensional (2D) version of one 
widely used in manganites: 
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Here Si is the localized spin operator representing the 
Ru +5 ion at site i, c ia is the operator creating an itinerant 



electron of spin a at this site, and Sj = X) Q /3 c \a a a-P c iP 
gives the spin of this electron. For simplicity we take the 
Hund coupling, Jh — * +oo. 

The model is solved using a classical Monte Carlo (MC) 
procedure for the localized spins, in conjunction with e: 
act diagonalization of the conduction electron system 
The localized spins are taken to be classical and of mag- 
nitude one. The conduction electrons are assumed to oc- 
cupy a single orbital, and from the condition Jh — > +oo 
only one spin projection is possible at each site. Then, 
for each configuration of localized spins, the effective hop- 
ping Uj of itinerant electrons between two sites i and j is 
affected by the factor (|i | Ti)i where | is the state of 
a spin 1/2 pointing in the direction of the localized spin 
j. Then: 
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where 0i and fa are the polar angles of spin Si. The 
resulting electronic energy levels are then filled by the 
available number of electrons in the canonical ensemble. 
The MC simulation proceeds from the partition function 
with classical spins. The general features of this model 
as a function of K/t and doping have bfifn studied be- 
fore in the context of the manganites. ffl The calcula- 
tions reported here are in a square of N =64 atoms with 
periodic boundary conditions and six mobile electrons 
(n = 6/64 = 0.1). We have taken K/t = 0.06. We ex- 
pect that the order of magnitude of t is near 0.25 eV, and 
then K ~ 180 K, near the observed Tm- 

We have calculated the spin structure factor defined as 
S(q) = J2i. 3 <S l .S j >e icl -^^yN 2 for all nonequivalent 
values of q in our 8x8 cell. In Fig. 1 we represent the 
two superlattice peaks of highest intensity as a function 
of temperature. Clearly, the peak at q = (tt, n) is the. 
dominant one, in agreement with neutron experiments.u 
An important fact is that at T = the amplitude of 
5(q) is smaller than one. This is consistent with experi- 
ments because the magnetic moment deduced from mag- 
netic measurementso is larger than that deduced from 
the amplitude of the superlattice reflections!] This reduc- 
tion in our simulations is due to the fact that in general, 
there are small ferromagnetic islands (magnetic polarons) 
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around the mobile electrons. Since these islands are not 
entirely compensated in our small systems and that their 
kinetics is slow at low temperatures, they give rise to a 
small contribution to S(0, 0) at small temperatures, as 
displayed in Fig. 1. Since we are dealing with a strictly 
2D system in our simulations (instead of 3D in the real 
system), and a continuous symmetry cannot be broken 
in 2D, Fig. 1 does not show a phase transition to the 
paramagnetic phase. 
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FIG. 1. Spin structure factor as a function of temperature 
for two wave vectors q: q = (tt, 7r) (diamonds), q = (0,0) 
(squares). 

The mobility of the ferromagnetic polarons makes it 
easy to orient them under small applied magnetic fields. 
Then, as observed experimentallyu the magnetization in- 
creases abruptly with small magnetic fields (~ 5 T) com- 
pared with the Neel temperature Tm — 133 K. This fact 
is unusual for an ordinary AF system. The magnetization 
curve of the localized spins in our system is shown in Fig. 
2. In spite of the difficulties with the sluggish statistics 
of the ferromagnetic polarons mentioned above, one can 
see an abrupt increase at low fields, due to the alignment 
of these polarons, followed by a linear increase with a 
smaller slope, in agreement with experiment.Q. 

Another unexpected result for an AF system is the 
fact that the inverse magnetic susceptibility at high 
temperatures extrapolates to a positive temperature for 
X~ x = 0. As shown in Fig. 3, our model also re- 
produces this result. The constant obtained fitting 
X = C/(T - 6) is 9 = (0.041 ± 0.004)t. This is a subtle 
issue which is related with the way in which the ferro- 
magnetic polarons are reduced as the temperature is in- 
creased, and their effect on the magnetic correlations on 
their neighborhood. A mean field picture in which the 
kinetic energy of the carriers is neglected seems unable to 
explain this result even qualitatively. The resulting value 
of the slope gives C = (0.314 ± 0.006)/z 2 , where \i is the 
magnetic moment of the localized spins. This value is 
near to the classical value /J? /3 for noninteracting spins. 
The value of C in the real system depends on the sain of 
the configurations 4d 5 and 4d 6 of Ru +5 and Ru+ 4 .□ 
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FIG. 2. Magnetization as a function of applied magnetic 
field at temperature T = 0.Q2K. M sat is the saturation mag- 
netization. 




FIG. 3. Inverse magnetic susceptibility as a function of 
temperature. The straight line is a fit using x^ 1 = (T— 0)/C. 

To conclude, we have proposed a double exchange 
model for RuSr2(Eu,Gd)Cu20s which is able to explain 
the main unusual features of these systems. Previous 
studies of the modcltltj indicate that for small densities 
there is phase separation in the model. However, in 
the real system, particularly due to the small number 
of carriersu, the effect of long-range Coulomb repulsions, 
not included in the model should be present and avoid 
macroscopic phase separation,^ while phase separation 
in a microscopic scale does not affect our conclusions. 
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